Doping dependance of the spin resonance peak in bilayer high- Tc superconductors 
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Motivated by a recent experiment on the bilayer Yi_a;Caa;Ba2Cu30j; superconductor and based 
on a bilayer t — J model, we calculate the spin susceptibility at different doping densities in the 
even and odd channels in a bilayer system. It is found that the intensity of the resonance peak in 
the even channel is much weaker than that in the odd one, with the resonance position being at 
a higher frequency. While this difference decreases as the doping increases, and both the position 
and amplitude of the resonance peaks in the two channels are very similar in the deeply overdoped 
sample. Moreover, the resonance frequency in the odd channel is found to be linear with the critical 
temperature Tc, while the resonance frequency increases as doping decreases in the even channel and 
tends to saturate at the underdoped sample. We elaborate the results based on the Fermi surface 
topology and the d-wave superconductivity. 
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I. INTRODUCTION 

Inelastic neutron scattering (INS) experiments have 
been playing an important role in the studies of the spin 
dynamics of high- Tc superconductors. They can provide 
direct information of the momentum and frequency de- 
pendence of the dynamical spin susceptibility. Over the 
past decade, one of the most striking features observed in 
the INS experiments is the resonant spin excitation. The 
resonance peak, which has been found in several classes of 
cuprate materials [U, S S 0, Q , has attracted much exper- 
imental and theoretical attention. This peak is centered 
at the momentum (tt, tt), with its intensity decreasing 
rapidly when the frequency moves away from (tt, tt). The 
resonance frequency is found to be in proportional to the 
critical temperature [1, @] ■ Theoretically, the origin of 
the spin resonance and its role on superconductivity are 
still open questions 0, i, i, [13, [H El, El ■ It has been 
proposed that the spin resonance is a collective spin exci- 
tation mode [H, El El, El E3, El El, 13 . Based on this 
scenario, many properties of spin fluctuations observed 
in the INS experiments have be explained consistently. 

In the YBazCuaOj, (YBCO) and Bi2Sr2CaCu208+:r 
(Bi-2212) family, one unit cell contains two Cu02 planes. 
The electronic states in different Cu02 layers belonging 
to one cell are strongly coupled at all doping levels. Thus, 
two modes of magnetic excitation are expected to exist in 
the bilayer systems, i.e., one in the even channel and the 
other in the odd channel according to the syrnmetry with 
respect to the exchange of the layers [111, [2I, [H, [31 • In 
earlier experiments, this expectation was only confirmed 
in the insulating YBCO samples [ll]. In the supercon- 
ducting state, the spin resonance mode was not observed 
in the even channel, presumably due to a much weaker 
intensity in this channel. Recently, the instrumentation 



advances have made it possible to resolve weaker fea- 
tures in the INS experiments. Two distinct resonance 
modes were observed in the superconducting state of bi- 
layer (Y,Ca)Ba2Cu30y samples [H, [l^. It was found 
that the resonance peak intensity in the even channel 
(J*^) is much weaker than that in the odd channel (/°), 
and the resonance frequency is higher than that of the 
odd channel. Very recently, the doping evolution of the 
resonance peak in both the even and odd channels of 
(Y,Ca)Ba2Cu30y was studied in detail by the INS ex- 
periments [ll]. In the overdoped samples {y — 7), the 
resonance position of the odd channel is close to that of 
the even channel. At this doping level, the two resonance 
modes have also closer intensities {1^/1° = 0.4). When 
the doping density decreases, the doping evolution of the 
resonance frequency in the odd channel seems to follow 
a similar doping dependence as Tc, while the resonance 
frequency seems to keep increasing in the even channel 
as the doping decreases and saturates to a constant in 
the underdoped sample. Moreover, the resonance peak 
intensity in the odd channel is also much larger than 
that in the even channel in the underdoped sample. The 
intensity ratio 1*^/1° decreases monotonously as doping 
decreases and reaches 0.05 in strongly underdoped sam- 
ples. 

Motivated by these experimental observations, we here 
present a detailed investigation of the doping dependence 
of the spin resonance mode in the even and odd channels. 
Following Ref. [1^, we employ a bilayer t — J type Hamil- 
tonian including the interlayer hopping and interlayer ex- 
change coupling. In order to examine the robustness of 
the doping dependance of the difference between the two 
resonance modes, we also look into in detail the effect of 
the interlayer hopping parameters t± and the interlayer 
exchange coupling J±. 

The article is organized as follows. In Sec. II, we 
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introduce the model and work out the formalism. In Sec. 
Ill, we perform numerical calculations and discuss the 
obtained results. Finally, we give a brief summary in 
Sec. IV. 



II. HAMILTONIAN AND FORMALISM 

We start with a Hamiltonian which describes a system 
with two layers per unit cell. 
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where I — 1,2 denotes the layer index. In the slave- 
boson approach, the physical electron operators cf^ are 
expressed by slave bosons 6-'^ carrying the charge and 
fermions representing the spin, cf^ 
mean-field level, we consider the order parameters A 
{.ff^ffl - /ff /]?) = ±Ao, (± depend on if the bond {ij) 



6f)t;^).Atthe 

(0 _ 



is in the x or the y direction), x^ij = J^aifia^ f- 
In the superconducting state, bosons condense b, 
h^P >= \/S, where(5 is the hole concentration. 
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Then, the mean-field Hamiltonian can be written as, 
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with 
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ASt' cos kx cos ky — /i, Ak = 2J'Ao(cosfca; — cosfcj^), 
£0 = 4iVJ'(xo + A§), and J' = 3 J/8. Diagonalizing the 
Hamiltonian, we can get the antibonding band (A) and 

bonding band (B) with the dispersion ^k^'^^ = ± t±. 
Here we use the momentum independent interlayer 
hopping constant tj_, being consistent with the recent 
angle resolved photocmission experiment on YBCO [30[, 
which reveals an obvious bilayer splitting along the 
nodal direction. 

The bare spin susceptibility can be expressed as. 
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FIG. 1: Panels (a) and (b) are the imaginary parts of the 
spin susceptibility versus the frequency for the wave vector 
Q = (7r,7r) with t± = O.IJ, J± = 0.15 J in the odd and even 
channels, respectively. Panel (c) is the resonance position as 
a function of the doping. Panel (d) is the intensity ratio of 
the spin susceptibilities between the even and odd channels. 
The quasiparticle damping F = 0.01. 



Here a, a' = A, B, E, 



(a) 



^^"^^ + A? is the quasipar- 



^k — V ''k "k 
tide energy, and /(w) is the Fermi distribution function. 

The bare even and odd channel spin susceptibilities 
which come respectively from the intraband and inter- 
band electronic transitions, are given by 



Xo(q,w) 



x(^-^)(q,c.)+x(^^^Hq,c.). (4) 



By including the correction of the antiferromagnetic 
(AF) spin fluctuations to the spin susceptibility in the 
form of the random-phase approximation (RPA), the 
renormalized spin susceptibilities for the even and odd 
channels can be obtained as 



X^(°)(q,c.) 



Xo (q,^) 



l + (5Jq± Ji) xxo'°^(q,c^)/2 



(5) 



where ± signs represent the spin susceptibility in the even 
and odd channels, respectively, J± is the interlayer ex- 
change coupling, Jq = J(cosga; -l-cosgj,) is the intralayer 
exchange. We here also include a to set the AF insta- 
bility at (5 = 0.02 [l3|. The mean- field order parameters 
Xo, Aq together with the chemical potential /i for dif- 
ferent doping 6 can be obtained from the self-consistent 
equations. The other parameters we choose are t = 2J, 
t' = -0.45<. 

Before we present our results, we wish to point out 
that the above formulas represent the spin susceptibility 
of the fermions. The spin susceptibility for physical elec- 
trons should be due to the boson condensation in the 
superconducting state. 
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FIG. 2: The resonance frequency as a function of the critical 
temperature Tc with tx_ — O.IJ and J± — 0.15J. 



III. RESULTS AND DISCUSSION 
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The odd and even channel spin susceptibilities at dif- 
ferent doping densities are shown in Figs. 1(a) and (b), 
respectively. As seen, the peak intensity in the odd chan- 
nel is sensitive to the doping, namely, it increases dra- 
matically as the doping decreases. While in the even 
channel it increases slowly as the doping decreases. To 
see the doping evolution more clearly, we plot the peak 
positions of the even and odd channels as a function of 
the doping in Fig. 1(c), and the ratio of the intensities 
of the spin resonance between the even and odd chan- 
nels, rj = 1^/1° versus doping in Fig. 1(d). From Fig. 1(c), 
the resonance frequency in the odd channel increases as 
the doping increases and saturates at the optimal dop- 
ing, then it decreases slightly in the overdoped regime. 
While in the even channel, the resonance frequency in- 
creases as the doping decreases, so that the peak posi- 
tions of the even and odd channels are closer and closer 
as the doping density increases, and the corresponding 
resonance frequencies are almost the same in the deeply 
overdoped sample. The difference of the intensities in 
the even and odd channels increases as the doping den- 
sity decreases, as seen in Fig. 1(d). The intensity ratio is 
only 0.1 in the strongly underdoped sample, and around 
0.4 in the overdoped region. Our results are qualitatively 
consistent with the experimental results [1^ . We also ex- 
amine the relationship between the resonance frequency 
and the critical temperature Tc by using an empirical for- 
mula Tc/Tcmax = 1 — c((5 — 0.16)^, whcrc c = 51 is used 
to ensure the AF instability to occur dX 5 = 0.02. The 
resonance frequency in the odd and even channels as a 
function of Tc is plotted in Fig. 2 [5 < 0.16). The reso- 
nance frequency (uj") is found to be proportional to Tc 
in the odd channel, which is in good agreement with the 
experimental results [ij, 5, 6]. While in the even channel, 
the spin resonance frequency (w^) depends weakly on Tc 
in the strongly underdoped sample, being also qualita- 
tively consistent with the very recent experiments [28| . 
in which the resonance frequency in the even channel is 
observed to increase in the overdoped sample and to sat- 



FIG. 3: (Color on line) (a) The intensity ratio of the spin 
resonance between the even and odd channels for different in- 
terlayer exchange coupling with t± — O.IJ. (b) The intensity 
ratio of the spin resonance between the even and odd channel 
for different interlayer hopping constants with J± — 0.15 J. 

urate to a constant in the underdoped sample. 

We now address the dependence of the intensity dif- 
ference of the spin resonance between the two channels 
on the parameters t±_ and J±. The intensity ratio as 
a function of the doping for different J± is plotted in 
Fig. 3(a). As seen, when J± increases, the intensity ratio 
decreases, indicating that the interlayer exchange cou- 
pling can strongly affect the even and odd channels and 
enhance the difference. However, the intensity ratio de- 
creases as the doping decreases for all J± we considered, 
indicating that our results presented above are robust 
against the variation of the interlayer exchange coupling 
J±. On the other hand, we can also see from Fig. 3 (a) 
that even if J_l = 0, the intensity in the odd channel is 
still significantly stronger than that in the even channel, 
suggesting that the interlayer exchange coupling is not 
the only contribution for the difference between the two 
channels. In fact, the other contribution comes from the 
interlayer single-particle hopping. To show this, we plot 
the ratio 77 versus the doping for different hopping con- 
stants t± in Fig. 3(b). As t± increases, the ratio decreases. 
So, the interlayer single particle hopping also contributes 
to enhance the difference between the two channels. Let 
us also consider the cases that t± depends on the doping 
density {t± cc S) and momentum [t± oc (cosfc^; — cosfcy)^]. 
Note that, the bilayer splitting is found to be momentum 
dependence in Bi-2212 bilayer systems [sT'l, without ob- 
serving the bilayer splitting in the nodal direction. We 
also examined that our results are robust for different 
types of the interlayer hopping constant t±, as shown in 
Fig.3(b). 

Now we elaborate the origin of the above features based 
on the Fermi surface topology. In Fig. 4 [s^l, we plot the 
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FIG. 4: (Color on line) The bare spin susceptibilities versus 
the frequency uj at diflterent doping densities in the odd and 
even channels with the quasiparticle damping F = 0.002. The 
dotted lines denote the real parts and the dashed lines imag- 
inary parts, respectively. The solid lines are the real parts of 
the RPA factor [1 + (qJq ± J_L)ReXo^°-'(Q, a;)/2] (scaledx2) 
with ± signs for the even and odd channels, respectively. 



imaginary and real parts of the bare spin susceptibility 
Xo^^ in the odd and even channels with different doping 
densities. We first address the spin excitation in the odd 
channel. As shown in Fig.4(a-c), the imaginary part of 
the bare spin susceptibility approaches to zero as the fre- 
quency is below the spin gap. At the edge of the spin gap, 
it has a step-like rise which arises from the flat band near 
(tt, 0) (Van Hove singularity). As a result, the real part of 
the bare spin susceptibility Re^o develops a sharp struc- 
ture near the edge. Consequently, a pole occurs when 
the real part of the RPA factor 1 -I- (5Jq — J_L)Rexg/2 
is equal to zero at the frequency u}° (slightly below the 
gap edge) and in the meantime the imaginary part of the 
spin susceptibility at uj° approaches to zero due to the 
spin gap. This suggests the formation of a spin collec- 
tive mode, which is ascribed to be the spin resonance. 
We can also see from Fig.4(a-c) that the frequency, at 
which a step-like rise occurs, increases as the doping de- 
creases starting from the optimal doping, but it decrease 
slightly in the overdoped regime. This explains the dop- 
ing dependence of the resonance frequency as presented 
in Fig. 1(c). On the other hand, the real part of the bare 
spin susceptibility also increases with the decrease of dop- 
ing. This leads to the pole position to be more and more 
below the spin gap edge (due to the finite damping T, 
the imaginary part of Xo is not zero slightly below the 



spin gap) and consequently to an increase in the renor- 
malized spin susceptibility. For the even channel case, 
an obvious difference is seen from Fig.4(d-f) that there 
are two step-like rises, instead of one in the odd channel. 
These two step-like rises come from the particle-hole scat- 
tering in the B and A bands, respectively, because the 
spin susceptibility in the even channel is contributed by 
the intra-band A ^ A and B ^ B scatterings as shown 
clearly in Eq.(4). The scattering within the B band leads 
to the step-like rise in Xo a lower frequency Q. b , while 
that in the A band leads to a higher frequency rise at Q.a- 
Since the rise is larger at O^^, the corresponding enhance- 
ment of the real part of Xo is larger there. Thus, the spin 
resonance peak occurs near VLa in this case. We note 
that, due to a smaller vertex olJq + J± in this channel, 
the pole condition 1 -I- (SJq -I- J±)KexQ/'2 = could not 
be satisfied at a large doping range. In the meantime, 
the corresponding imaginary part of Xo is of appreciable 
value because of the scattering in the B band. So, the 
resonance peak in this case is basically a quasi-resonance 
peak, with its intensity being much lower than that in 
the odd channel as shown in Fig. 1(b). Also in contrast 
to the case of the odd channel, the frequency, at which 
the high-energy step-like rise occurs, increases with the 
decrease of doping. Therefore, the spin resonance fre- 
quency at the even channel increases upon reducing the 
doping as shown in Fig. 1(c). On the other hand, one 
can see from Fig.4(d-f) that both the real and imaginary 
parts of Xo i^ the lower side of CIa do not change much 
with doping. Considering the appreciable increase of the 
resonance intensity with the decrease of doping in the 
odd channel, it is expected that the difference between 
the spin resonance peak intensities of the two channels 
increases as the doping decreases. 

These features can be traced to the evolution of the 
Fermi surface with doping. We present the normal state 
Fermi surface in Fig. 5. As discussed above, the step-like 
rise is near the spin gap edge. In the zero temperature 
limit, the bare spin susceptibility [Eq.(3)] can be rewrit- 
ten as Imx("'"'Hq,w) oc - q))- Here, 
flica') (^^^ ^ ^u^+'^'^+q denotes the energy to break a 
pair and excite two quasiparticles from the superconduct- 
ing condensed state, and has a minimum of the exciting 
energy (MINk[ri'-"'" •* (k, q)]) when the wave vector q is at 
Q = (tt, tt) where the spin resonance is observed, which is 
just the spin gap. Because of the d-wave symmetry of the 
superconducting gap and energy band structure, the ex- 
citation within the A band with the minimum excitation 
energy is the Af'-to-M excitation as shown in Fig. 5 (a), 
while that within the B band and that of the interband 
correspond respectively to the N'-to-N and O'-to-O exci- 
tations, where M, M', N, N' are the crossing points (hot 
spot) of the Fermi surface with the magnetic Brillouin 
zone boundary, and 0{0') are the crossing points of the 
B{A) band Fermi surface with the (tt, tt) shifted images 
of the A{B) band Fermi surface. From Fig. 5(b), we can 
see that the hot spot of the B band moves towards the 
nodal direction as the doping decreases, and consequently 



5 




k 1% 

X 

FIG. 5; (Color on line) (a) The normal state Fermi surface 
at the doping 5 = 0.16 with t±_ = O.IJ. The bold solid and 
dashed lines are the Fermi surface of the A and B bands, re- 
spectively. The dotted and the dash-dotted lines in the first 
quadrant are the (tt, tt) shifted images of the A and B band 
Fermi surfaces in the third quadrant. The solid and dashed 
arrows denote the intraband and interband scatterings, re- 
spectively, (b) The Fermi surfaces of the A (solid lines) and 
B bands (dashed lines) in the first quadrant of the Brillouin 
zone at different doping densities. 

the magnitude of the corresponding superconducting gap 
decreases. The A band depends weakly on the doping 
density, but the magnitude of the superconducting gap 



increases with doping, as calculated from the mean-field 
theory here and observed in experiments [1^. Because 
the spin excitation in the odd channel comes from the 
O'-to-O excitation, its spin gap decreases with the de- 
crease of doping. While, the high-energy step-like rise in 
the even channel is contributed by the excitation from 
M' to M, it increases with the decrease of doping. 



IV. SUMMARY 

In summary, we have examined the doping evolution 
of the spin susceptibility in the even and odd channels in 
the bilayer high- Tc superconducting materials based on 
the bilayer t — J type model. In the bonding and anti- 
bonding band representation, there exist two channels of 
the spin excitation according to the intraband scattering 
and interband scattering. Each channel has its distinct 
resonant mode. In the odd channel, i.e., the interband 
scattering, the spin susceptibility shows a strong dop- 
ing dependence. As the doping decreases, the intensity 
increases dramatically and the resonance frequency is lin- 
ear with Tc- The resonance frequency in the even chan- 
nel approaches to that in the odd channel and the ratio 
between the two channels is around 0.4 in the overdoped 
region. As the doping decreases, the resonance frequency 
increases and saturates at the strongly underdoped sam- 
ple. In addition, it has been found that the differences of 
the resonance positions and intensities between the two 
channels are enlarged as the doping decreases. Our re- 
sults are well consistent with the experiments. We have 
elaborated the results based on the topology of the Fermi 
surface and the d-wave superconductivity. 
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